Ultra-High Energy Cosmic Ray Probes of Large Scale Structure and Magnetic Fields* 



Giinter SigP^'', Francesco Miniati'^, Torsten A. Enfilin'^ 
" GReCO, Institut d'Astrophysique de Paris, C.N.R.S., 98 bis boulevard Arago, F-75014 Paris, France 
Federation de Recherche Astroparticule et Cosmologie, 
Universite Paris 1, 2 place Jussieu, 75251 Paris Cedex 05, France 
Max-Planck Institut fiir Astrophysik, Karl-Schwarzs child- Str. 1, 85741 Garching, Germany 

We study signatures of a structured universe in the multi-pole moments, auto-correlation function, 
and cluster statistics of ultra-high energy cosmic rays above 10^^ eV. We compare scenarios where 
the sources are distributed homogeneously or according to the baryon density distribution obtained 
from a cosmological large scale structure simulation. The influence of extragalactic magnetic fields is 
studied by comparing the case of negligible fields with fields expected to be produced along large scale 
shocks with a maximal strength consistent with observations. We confirm that strongly magnetized 
observers would predict considerable anisotropy on large scales, which is already in confiict with 
current data. In the best fit scenario only the sources are strongly magnetized, although defiection 
can still be considerable, of order 20° up to 10^" eV, and a pronounced GZK cutoff is predicted. 
We then discuss signatures for future large scale full-sky detectors such as the Pierre Auger and 
EUSO projects. Auto-correlations are sensitive to the source density only if magnetic fields do 
not significantly affect propagation. In contrast, for a weakly magnetized observer, degree scale 
auto-correlations below a certain level indicate magnetized discrete sources. It may be difficult 
even for next generation experiments to distinguish between structured and unstructured source 
distributions. 

PACS numbers: 98.70.Sa, 13.85.Tp, 98.65.Dx, 98.54.Cm 



I. INTRODUCTION 

The origin of ultra-high energy cosmic rays (UHECRs) 
is still one of the most challenging problems of mod- 
ern astrophysics. It is an open question which mecha- 
nism is responsible for producing particles up to 10^" eV 
and beyond and where the corresponding sources can be 
found [ll |2| ■ Although statistically meaningful informa- 
tion about the UHECR energy spectrum and arrival di- 
rection distribution has been accumulated, no conclusive 
picture for the nature and distribution of the sources 
emerges naturally from the data (for a short overview 
on the relevant literature see Ref. ^J). There is on 
the one hand the approximate isotropic arrival direc- 
tion distribution which indicates that we are observ- 
ing a large number of weak or distant sources. On the 
other hand, there are also indications which point more 
towards a small number of local and therefore bright 
sources: First, there seem to be statistically significant 
multi-plets of events from the same directions within a 
few degrees 0, 0. Second, nucleons above ~ 70EeV 
suffer heavy energy losses due to photo-pion production 
on the cosmic microwave background — the Greisen- 
Zatsepin-Kuzmin (GZK) effect p — which limits the 
distance to possible sources to less than ~ 100 Mpc 0- 
For a uniform source distribution this would predict a 
"GZK cutoff" , a drop in the spectrum. However, the 
existence of this "cutoff" is not established yet from the 
observations In fact, whereas a cut-off seems con- 
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sistent with the few events above 10^" eV recorded by 
the fluorescence detector HiRes 9], it is not compatible 
with the 11 events (also above 10^*^ eV) measured by the 
AGASA ground array ^1^]. The solution of this problem 
may have to await the completion of the Pierre Auger 
project jllj which will combine the two complementary 
detection techniques adopted by the aforementioned ex- 
periments. 

Such apparently contradicting hints could easily be 
solved if it would be possible to follow the UHECR tra- 
jectories backwards to their sources. However, this may 
be complicated by the possible presence of extragalactic 
magnetic fields, which would deflect the particles during 
their travel. Furthermore, since the GZK-energy losses 
are of stochastic nature, even a detailed knowledge of the 
extragalactic magnetic fields would not necessarily allow 
to follow a UHECR trajectory backwards to its source 
since the energy and therefor the Larmor radius of the 
particles have changed in an unknown way. Therefore it 
is not clear if charged particle astronomy with UHECRs 
is possible in principle or not. And even if possible, it 
remains unclear to which degree the angular resolution 
would be limited by magnetic deflection. 

Quite a few simulations of the effect of extragalactic 
magnetic fields (EGMF) on UHECRs exist in the htera- 
ture, but usually idealizing assumptions concerning prop- 
erties and distributions of sources or EGMF or both are 
made: In Refs. T?, T?, TI, T?, T?] sources and EGMF fol- 
low a pancake profile mimicking the local supergalactic 
plane. In Ref. |l3| highly structured EGMF have been 
obtained from constrained simulations, but the source 
distribution has been assumed homogeneous. In other 
studies EGMF have been approximated in a number of 
fashions: as neglig ible 111 ll9j, as uniform '2(|, or as 
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organized in spatial cells with a given coherence length 
and a strength depending as a power law on the local 
density 

However, the presence of the above mentioned appar- 
ently contradicting hints indicate that the existing data 
set might already carry information on non-trivial prop- 
erties of sources and EGMF. Here, we want to address 
the following questions relevant to charged particle as- 
tronomy: 

1. Do we observe a large number of dim or a small 
number of bright UHECR sources? 

2. Is the source location distribution statistically ho- 
mogeneous, or does it follow the matter distribution 
in the local Universe? 

3. Are the particles strongly deflected by intergalactic 
magnetic fields? 

4. What are the magnetic fields surrounding sources 
and observer ? 

5. Can we discriminate between the case of sources 
with practical identical luminosities and the case 
of a power-law luminosity distribution? 

The means by which we want to study these questions 
are comparisons of simulated datasets to the observed 
one by statistical tests on the arrival direction multi-pole 
moment and auto-correlation distributions, on the multi- 
plet statistics, and on the UHECR energy spectrum. The 
simulated UHECR events are produced by following the 
trajectories of particles through a large-scale structure 
simulation which included a numerical model for the gen- 
eration and evolution of the EGMF. We thereby extend 
our former studies .'i. 22] to a larger parameter space and 
a higher degree of realism. 

The most important results of these former studies are: 
A local component of sources within ^ 100 Mpc alone 
cannot explain satisfactorily the observed isotropy at en- 
ergies < 4 X lO^^eV |3l. In Ref. [13 and the present 
study we therefore take into account sources at cosmo- 
logical distances by periodically repeating the large scale 
structure simulation box. As a result we found that in 
combination, a comparison of spherical multi-poles for 
I < 10 and of the auto-correlation at angles 9 < 20° 
between observed and simulated data moderately favors 
a scenario in which (i) UHECR sources have a density 
Us ~ 10~^ Mpc"'^ and follow the matter distribution (ii) 
magnetic fields are relatively pervasive within the large 
scale structure, including filaments, and with a strength 
of order of a fiG in galaxy clusters (iii) the local extra- 
galactic environment is characterized by a weak mag- 
netic field below 0.1 /iG. This is in contrast to Ref. j3| 
where the neglect of cosmological sources marginally fa- 
vored observers immersed in ~ 0.1 fj,G fields. Finally, 
we found that the degree-scale auto-correlation functions 
above ~ 4x lO^^eV can serve as a discriminator between 
magnetized and unmagnetized sources. 



In the present paper we specifically focus on signatures 
for the structure, density and luminosity distribution of 
UHECR sources, as well as for magnetic fields surround- 
ing sources and observer. We will find, as expected, that 
for most of these observables, the current data set does 
not allow to clearly distinguish between limiting cases. 
We will therefore discuss how future experiments such as 
the Pierre Auger ^nd EUSO [HI projects will im- 
prove the prospects to measure these observables. 

Here, we restrict ourselves to UHECR nucleons, and 
we neglect the Galactic contribution to the deflection 
of UHECR since typical proton deflection angles in 
galactic ma gne tic fields of several /uG are < 10° above 
4 X 10^^ eV l2J|, and thus in general are small compared 
to extragalactic deflection in the scenarios studied in the 
present paper. 

The simulations are described in more detail in the 
next section. There we also describe the general fea- 
tures of our method and define the statistical quantities 
used for comparison with the data. Sect. 3 presents how 
large scale multi-poles and small-scale auto-correlations 
probe magnetization and UHECR source characteristics. 
It constitutes the main part of the present paper. We 
conclude in Sect. 4. 



II. OUTLINE OF THE NUMERICAL 
TECHNIQUE 

In the next section we will investigate, from a sta- 
tistical point of view as defined below, the viability of 
various scenarios for the propagation of UHECRs in an 
extragalactic environment. These scenarios are listed in 
Tab.^and differ in the UHECR source distributions, the 
strength of the EGMF, and the location of an hypotheti- 
cal observer on Earth. The latter will also lead to differ- 
ent strengths of the EGMF within a few Mpc from the 
observer. In the following sections we describe in detail 
how the various scenarios are characterized. 



A. Large Scale Structure and Extragalactic 
Magnetic Fields 

The magnetized extragalactic environment which we 
use for our experiments is produced by a simulation of the 
large scale structure of the Universe. The simulation was 
carried out within a computational box of 50ft.^^Mpc 
length on a side, with normalized Hubble constant h = 
iJo/(100 km s~^ Mpc~^) = 0.67, and using a comoving 
grid of 512'^ zones and 256^ dark matter particles. This 
is the same large scale structure simulation that was used 
in Refs. ^^'^ further described in 25]. 

The EGMF adopted here is based on the numerical 
model developed in [2^ which has been shown to be 
compatible with existing Faraday rotation measures with 
lines of sight both through clusters and the diffuse inter- 
galactic medium. Thus, at simulation start the EGMF 
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FIG. 1: Distribution of the number of cells n as a function of 
the maximum M of the Alfven Mach number and the plasma 
beta parameter. 



was initialized to zero and subsequently its seeds were 
generated at cosmic shocks through the Biermann bat- 
tery mechanism \2T\ . This approach is alternative to the 
case in which the initial magnetic field is set uniform 
over the whole simulated volume. Since cosmic shocks 
form primarily around collapsing structures including fil- 
aments, the above approach avoids generating EGMF 
in cosmic voids. An alternative, more realistic but also 
much more complicated scenario is discussed in Ref . |28| , 
in which magnetic fields are injected into the intergalac- 
tic medium by galactic outflows. Whichever the mecha- 
nism that generates it, the magnetic field is then evolved 
according to the induction equation and is therefore am- 
plified in different parts of the universe by shear flows 
and compression according to the velocity field provided 
by the simulated gas component. 

As already pointed out in Ref. 0, given the tiny val- 
ues of the initial seeds generated through the Biermann 
mechanism, the field strength at simulation end (when 
the cosmological redshift equals zero) is much smaller 
than what is observed in galaxy clusters via a number of 
experiments [2^. Therefore, in order to avoid such dis- 
crepancies, it is necessary to change the normalization of 
the simulated magnetic field strength. Our renormaliza- 
tion procedure simply involves a rescaling of the overall 
magnetic field in the computational box, such that the 
magnetic field in the core region of a Coma-like galaxy 
cluster is predicted to be of order of a uG or so, as in- 
dicated by Faraday rotation measures |23|. As a result 
of this rescaling, the magnetic field strength volume av- 
eraged over ~ 0.5 Mpc within typical cluster cores is be- 
tween 0.7 and 2.5 fxG. Since the magnetic field strength. 
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FIG. 2: The cumulative filling factors for EGMF strength, 
(middle panel) baryon density (lower panel, in units of av- 
erage baryon density) above (decreasing curves) and below 
(increasing curves) a given threshold, as a function of that 
threshold. 

B, in collapsed structures follows a well defined scaling 
relation with the structures virial temperature [30j |. the 
renormalization of the magnetic field as described above 
can be easily carried out even though Coma like galaxy 
clusters do not form in our simulation due to the rela- 
tively small computational box. 

Lacking direct measurements of magnetic fields in fila- 
ments, we assume that the topology and relative strength 
of the large scale intergalactic magnetic field in different 
parts of the Universe is as reproduced by our numerical 
simulation. However, as already pointed out the result- 
ing EGMF is consistent with statistics of existing Fara- 
day rotation measures with line of sight through filaments 
despite the fact that the magnetic fleld strenght can be 
close to the equipartition value with the gas total energy 

m. 

It is worth pointing out that, because of the very scarce 
observational constraints on intergalactic magnetic flelds, 
both our assumptions have limitations. For example, an 
independent experiment based on the detection from the 
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magnetic field 




FIG. 3: Log-scale two-dimensional cut through magnetic field 
total strength in Gauss (color scale in Gauss, upper panel) and 
baryon density in units of average baryon density (color scale, 
lower panel). The observer is in the center of the figures and 
is marked by a star. The EGMF strength at the observer is 
~ 0.1 jj,G. Note that both panels correspond to the same cuts 
through the full large scale simulation box. 



Coma cluster of radio synchrotron emission and hard X- 
rays interpreted as inverse Compton emission, would sug- 
gest that magnetic fields in this cluster are lower by about 
an order of magnitude compared to what we are assum- 
ing [sij . Similarly, and more importantly, one cannot 
exclude the possibility of much smaller magnetic fields 
in filaments than we assume, although some evidence for 
magnetic fields in filaments at the level of a tenth of a 
fiG may already exist 32] . 

According to our simulation scenario and with the 
assumptions made, we find that EGMF are significant 
only within filaments and groups/clusters of galaxies. In 
Fig. n we test the assumption that the magnetic field 
is passive. In order for the magnetic field to be dy- 
namically important both the Alfvenic Mach number. 
Ma — v/vA — f /-B/(47rp)^/^, and the plasma (3 parame- 
ter, (3 — Pgas/ {B'^ /8tt) must be smaller than unity. Thus 
in Fig. n we plot a histogram of the fraction of cells as a 
function of M = max(M^, /?). The histogram shows that 
the condition of dynamically unimportant magnetic field 
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FIG. 4: Same as Fig.|21 but for an observer situated in a small 
void where the EGMF strength is ~ 10~" G. 



is violated only in a very small fraction of the volume, 
which does not affect the evolution of the simulation in 
any significant way. Furthermore, only a fraction of the 
cells in those bins is characterized by a a magnetic field 
capable of affecting the trajectory of UHECRs. In any 
case, as discussed in Sec. IIVI we consider the case of a 
field strength normalization reduded by a factor 10 with 
respect to our fiducial model. This would correspond to 
a left-shift of the x-axis in Fig. of two decades, such 
that the EGMF is virtually dynamically unimportant in 
all cells of the simulation. 

About 90 percent of the volume is filled with fields 
< 10 nG and in the voids fields are < 10"" G. In Fig.H 
we show the cumulative filling factors for both the EGMF 
strength and the baryonic gas. The volume field averages 

are {\B\) ~ 1.5 x 10"^ G, and ((S^))^^^ ~ 7.9 x IQ-^G, 
with coherence lengths ^ Mpc in the strong field regions, 
and thus compatible with Faraday rotation bounds 2^ 
|2^ . Note that the hallmark of a highly structured field 

1/2 

is a ratio ((S^)) / {\B\) :$> 1, as in the present case. 

As in Refs. [1,|23|, we explore the case of two observers 
located at different positions within the simulated vol- 
ume: The first is located in a filament-like structure with 
EGMF ^ 0.1 /iG, and the second at the border of a small 
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void with EGMF - 10"" G. Figs. H and H show two- 
dimensional cuts through the EGMF and baryon den- 
sity around these two observer positions. Notice that the 
structure of the magnetic field is quite more extended 
than that for the baryonic gas. Thus a UHECR pro- 
duced by a source where matter density is high will be 
subject to the action of magnetic fields within an ex- 
tended volume surrounding the source, before breaking 
into a void where magnetic fields are much weaker. A 
relatively large structure about 17 Mpc away from the 
weak field observer is identified for calculation purposes 
as the Virgo cluster. We orient our terrestrial coordinate 
system so that this cluster is close to the equatorial plane. 

B. UHECR Sources 

For a given number density of UHECR sources, Ug, we 
explore both the case in which their spatial distribution 
is either proportional to the local baryon density, as in 
Ref. or completely homogeneous. Further, due to the 
unknown source positions and properties, there will be a 
cosmic variance in the results. In order to be conservative 
and maximize this variance, we will therefore assume that 
all UHECR sources, are distributed in luminosity, Qi, so 
that their contribution of UHECR per log Qi is roughly 
constant with Qi. In addition we allow the spectral index 
tti of the emitted power-law distributions of UHECRs to 
vary and assume that each source accelerates UHECRs 
up to 10^^ eV. Specifically, our assumption can be sum- 
marized as follows, 

^ (X Qr2-2 foj. 1 < < Q,„^^ 
uL^i 

— const for a — Aa < ai < a + Aa (1) 

dat 

i^max - lO'^eV, 

where Qi has been put in dimcnsionlcss units. The power 
law distribution in Qi could be further motivated by the 
luminosity function of the EGRET 7— ray blazars which 
has this sha pe in the power range 10''^ erg s""'^ ^ Qi ^ 
lO^'^ergs"^ Constant source characteristics corre- 

spond to Qmax = 1, Aa = in Eq. If not exphcitly 
stated that sources with identical properties are assumed, 
we will always use Qmax = 100, Aa = 0.1. Finally, 
the actual value of a, representing the central value of 
the power law index distribution, together with the total 
power of injected UHECRs are left as free parameters to 
be obtained from a best fit analysis when the simulation 
results are compared with observational data. 

We also assume that neither total power, Qi, nor the 
power law spectral index, ai, change significantly on the 
time scale of UHECR propagation. For energies higher 
than the GZK cutoff this is not an issue because the 
sources are nearby and the propagation time, given by 
the distance in light years added by the time delay pro- 
vided in Fig. [SI is less than the typical duty cycle of, say. 
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FIG. 5: The distribution of UHECR detection energy E and 
delay time r with respect to straight line propagation, for 
scenario 6 in Tab.|l] averaged over 26 realizations of 10'* sim- 
ulated trajectories above lO'^^ eV each. 



a radio source. On the other hand, for UHECRs of lower 
energies the propagation time can be up to a few Gyr, see 
Fig. 13 However, the flux of UHECRs at these energies is 
dominated by many sources at relatively large distances 
and, therefore, should not be sensitive to time variations 
of individual sources. 

Taking advantage of the periodic boundary conditions, 
we build a "hyper volume" by adding periodic images of 
the simulation box until the linear size of the enclosed 
volume is 1.5h~^ Gpc, that is larger than the energy 
loss length of nucleons above lO'^^ eV. The distributions 
of sources and the EGMF are identical in each box for a 
given realization. Further, in each replicated box within 
the central Gpc of the hyper volume we place an observer 
at the same position (relative to the box). This ensures 
that each observer is still surrounded by several 100 Mpc 
of sources in each direction. Since the e nerg y loss length 
for nucleons is < 1 Gpc above 10^^ eV [sj, all relevant 
distant sources are taken into account in this manner. 



C. UHECR Propagation and Event Detection 

For each configuration, consisting of a choice of source 
positions and a set of parameters Qi and ai, many nu- 
cleon trajectories originating from sources in the Gpc^ 
volume were computed numerically. As in previous work 

Jiiiii J£jj particles are propagated taking into ac- 
count Lorentz forces due to the EGMF and energy losses. 
In this respect, pion production is treated stochastically, 
whereas pair production is treated as a continuous en- 
ergy loss process. Cosmological redshift cannot be taken 
into account because the propagation time is not known 
before hand. However, this is a minor effect since the 
travel time of the trajectories does not exceed ^ 3 Gyr, 
corresponding to z ~ 0.6 down to 10^^ eV, see Fig. [S] 
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Trajectories connecting sources and observers in differ- 
ent copies of the simulation box are taken into account. 

An event was registered and its arrival direction and 
energy recorded each time the trajectory of the propagat- 
ing particle crossed a sphere of radius 1.5 Mpc around 
one of the observers. Events at replicated observers in 
different boxes are recorded by the same counter. The 
counter is stopped when 10** events (a realization) have 
been recorded. For more details on this method see 
Refs. [Ullllll. Modeling the observer as a sphere 
corresponds to an average over observers located on 
that sphere. The concrete directions of features in the 
UHECR sky distribution can never be accurately recon- 
structed in this approach. However, quantities which 
only depend on angular distances such as the multi-poles, 
the auto-correlation function, and the number of event 
clusters can properly be accounted for by this approach. 

Having provided all the relevant definitions, for vari- 
ous scenarios explored in the next section we build sta- 
tistical samples in the following way. We generate ten 
realizations for the source positions. For each one of 
them we use the distributions in Eq. with Qmax — 
100, Aa = 0.1, and typically generate 50 realizations 
with different power Qi and injection spectral index a^, 
whereas the maximal acceleration energy is held constant 
at 10^^ eV for simplicity. Thus, for each scenario typically 
500 realizations corresponding to different configurations 
of source positions and emission characteristics are simu- 
lated. The variation of results with these configurations 
constitutes cosmic variance. 



D. Data Processing 

The process of converting the simulated events to 
quantities that can be compared with observations was 
described at length in Ref. j^]- Here we simply summa- 
rize the procedure. 

For each realization discussed above, the events were 
used to construct arrival direction probability distribu- 
tions above a given energy on regular grids of solid angle 
bins of size in both declination and right ascension. 
For the full sky, for example, there are ~ 180°/A6' bins 
in declination and, at declination 5, ~ (360°/A6') cos (5 
bins in right ascension. We chose A0 small compared to 
the instrumental angular resolution; typically A0 — 0.5°. 
The sky distributions are multiplied with the solid-angle 
dependent exposure function for the respective experi- 
ment and are convolved with the angular resolution which 
is ~ 1.6° above 4 x lO^^ eV and ~ 2.5° above 10^^ eV for 
the AGASA experiment 0. Due to the poorer angu- 
lar resolution of the SUGAR data [s^, we will use a 
resolution of 2.5° when combining AGASA and SUGAR 
data. Energy resolution effects, supposed to be of order 
AE/E ~ 30%, are also taken into account. 

For the exposure function U!{5) we use the parameter- 
ization of Ref. IT^ which depends only on declination 



oj{5) (X cos(ao) cos((5) sin(a. 



where a. 



with 



- am sin(ao) sin((5) , 

if ^ > 1 

if^<-l , (2) 
cos ^(^) otherwise 

cos(0„,) — sin(ao) sin((5) 
cos(ao) cos((5) 



For the AGASA experiment 10\ gq = 35°, 0^ = 60°, 
and for the former SUGAR experiment on the south- 
ern hemisphere cq = —30.5°, 9m — 55°. For a full-sky 
Pierre Auger type experiment we add the exposures for 
the Southern Auger site with aq = —35° and a putative 
similar Northern site with oq = 39°, and 9m — 60° in 
both cases, with an assumed angular resolution of ~ 1°. 

From the distributions obtained in this way typically 
200 mock data sets consisting of Nohs observed events 
were selected randomly. For each such mock data set or 
for the real data set we then obtained estimators for the 
spherical harmonic coefficients C{1), the auto-correlation 
function N{9), and the number of multi-plets M(n) of 
n events within an angle 9m- As in Refs. 0, ll^ the 
estimator for C{1) is defined as 



Oil) 



21 



m——l \ 2—1 / 



(3) 



where uji is the total experimental exposure at arrival 
direction , M — X^i^i" "^^^ snra of the weights 

l/cji, and Yimiu^) is the real- valued spherical harmonics 
function taken at direction . Also as in Refs. ^ 22^ the 
estimator for N{9) is defined as 



N{9) 



C 

W) 



E 



1 if 9ij is in same bin as i , 
otherwise ' ' ^ ' 



and S{9) is the solid angle size of the corresponding bin. 
In Eq. the normalization factor C — fle/{Nohs{Nohs~ 
1)), with ile denoting the solid angle of the sky region 
where the experiment has non- vanishing exposure, is cho- 
sen such that an isotropic distribution corresponds to 
N{9) = 1. 

Finally, the multi-plets are obtained as follows: For a 
given set of A'obs ordered events, the first event is con- 
sidered as a singlet. Then for all other events we check 
if they lie within distance 9m of one of the other events 
defining the center of a multi-plet. If yes, the multiplicity 
of that multi-plet is increased by one, if not, the event 
is defined as the center of a new multi-plet, starting as a 
singlet. 

In passing we note that other statistical quantities 
have been considered in the literature, such as cumu- 
lative distributions in right ascension and declination in 
the context of the Kolmogorov-Smirnov test , and the 
"information-dimension" of the sky distribution |3^ . We 
will not use these statistics in the present work. 
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The different mock data sets in tfie various realiza- 
tions yield the statistical distributions of C{1), N{9), and 
M{n). One defines the average over all mock data sets 
and realizations as well as two errors. The smaller er- 
ror (shown to the left of the average in the figures be- 
low) is the statistical error, i.e. the fluctuations due to 
the finite number iVobs of observed events, averaged over 
all realizations. The larger error (shown to the right of 
the average in the figures below) is the "total error" , i.e. 
the statistical error plus the cosmic variance. Thus, the 
latter includes the fluctuations due to finite number of 
events and the variation between different realizations of 
observer and source positions. 

Given a set of observed and simulated events, after 
extracting some useful statistical quantities S'i, namely 
C;, N{6), and M{n) defined above, we define 

Ef data '^'i.siniu \ /^\ 
I j 

where the index i runs over multi-pole angular bin of 
6, and multiplicity n, respectively. Here, refers to 

Si obtained from either the real data set or the simulated 
mock data sets, and S'i^simu and AS'i^simu are the average 
and standard deviations of the simulated data sets. Thus, 
there is a Xn for the real data, and a Xn for each of the 
simulated mock data sets which consist of the same total 
number of UHECR events. This measure of deviation 
from the average prediction can then be used to obtain an 
overall likelihood for the consistency of a given theoretical 
model with an observed data set by counting the fraction 
of simulated data sets with Xn larger than the one for the 
real data. 



III. RESULTS: PROBING EGMF AND UHECR 
SOURCE CHARACTERISTICS 

We now turn to a systematic discussion of signatures of 
magnetization and UHECR source characteristics in the 
angular power spectrum, the auto-correlation function, 
and the clustering of the UHECR arrival distributions. 

The scenarios studied are presented in Tab. ^together 
with a statistical measure of their likelihoods. To sum- 
marize, UHECR sources whose number density is given 
in column 2 are distributed either proportionally to the 
simulated baryonic density or homogeneously ("yes" or 
"no" , respectively, in column 3). The observer is either in 
a region with or without appreciable magnetic fields, as 
quantified more precisely in Fig.|21and^ respectively. Fi- 
nally, the EGMF is either taken from the simulation with 
local value as indicated in column 5, or completely ne- 
glected ( "no EGMF" ) . For all scenarios our statistical as- 
sessments are based on typically 200 simulated mock data 
sets each for 500 different realizations of source locations 
and emission characteristics, as discussed in Sects. H.C 
and H.D above. Scenarios 1-6 in Tab. ^ have already 
been presented in Ref. !2^, whereas scenarios 7-10 are 
new. 



Columns 7-14 of Tab. U show the likelihood signif- 
icances discussed in the previous section obtained by 
comparing the predictions from each scenario with cur- 
rently available experimental data. Particularly, as in 
previous studies, we carry out a comparison in terms 
of multi-poles, auto-correlations and multi-plet statis- 
tics. As for the experimental data we use the 57 
AGASA-I-Akeno events above 4 x 10^^ eV, when study- 
ing the auto-correlation function and clustering proper- 
ties of UHECRs which are sensitive to small scales. For 
the large scale multi-poles, however, larger sky coverage 
including the Southern hemisphere is desirable in order 
to get realistic estimates of the true multi-poles. There- 
fore, following Ref. 38], when comparing with multi- 
poles Z < 10 we will also use the combination of 50 
events observed by AGASA (excluding 7 events observed 
by Akeno) and 49 events seen by SUGAR for a total of 99 
events above 4x 10^^ eV. In fact, as pointed out by the au- 
thors of Ref. [Hi, the AGASA and SUGAR experiments 
have comparable exposure in the northern and southern 
hemisphere, respectively. In addition, while SUGAR'S 
angular resolution is much worse than for AGASA and 
in general prevents a combination of the two data sets, 
this does not affect multi-poles / < 10 because they are 
not sensitive to scales < 10°. Comparing columns 7 with 
10 and 8 with 11 in Tab. ^ shows that the inclusion of 
SUGAR data does not significantly change the likelihood 
ranking of simulated scenarios in terms of multi-poles. 

Finally, using the combined exposure of 
AGASA-I-SUGAR, for each scenario we compute 
the expected values of the multi-pole coefficient for the 
~ 1500 events observed above 10^^ eV. This is interesting 
because no sign of anisotropy was found by either the 
AGASA or the SUGAR experiment at these energies. 
Since data down to 10^^ eV are not publicly available, we 
simply compare our results with a completely isotropic 
distribution. The corresponding likelihoods are also 
summarized in Tab.^ 

As discussed in Ref. the current best fit is pro- 
vided by scenario 6, i.e. for structured sources of density 
Us ~ 10~^ Mpc~^ roughly following the baryon density 
and immersed in fields up to a few micro Gauss, whereas 
the observer is surrounded by fields ^ 10~^ G. In ad- 
dition, we note that if the sources are homogeneously 
distributed, the best fit density of the three values con- 
sidered in Tab. m (scenarios 5, 7, and 8) turns out to be the 
intermediate one, ^ 2.4 x 10~^, in rough agreement 
with Ref. However, the likelihood is a very shallow 
function of and implies uncertainties of at least an 
order of magnitude. 

To illustrate the general impact of an EGMF on propa- 
gation, we show in Fig.[Slthe distributions of delay times 
T with respect to straight line propagation with arrival 
energies i?, averaged over all realizations for scenario 6 
in Tab. Note that the various peaks at small time de- 
lays T < 10** yr are due to UHECRs from discrete, nearby 
sources which mostly see the relatively weak fields around 
Earth in this scenario. One can show that the number 
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TABLE I: List of simulated scenarios. The columns contain the number assigned to the scenario, the source density, whether 
the sources are distributed as the baryon density in the large scale structure simulation box or homogeneously (yes/no), the 
observer position 1 or 2 corresponding to Fig. |3and|ll respectively, the magnetic field strength at the observer location (zero 
indicates no fields, whereas a number indicates the EGMF obtained from the large scale structure simulation), the best fit power 
law index in the injection spectrum E~°', and the overall likelihoods of fits to the data. The first six likelihoods are for the 
multi-poles Eq. ^ above the energy indicated as superscript in EeV and over the range of I indicated as subscript. "AGASA 
only" and "AGASA+SUGAR" indicates which exposure functions and data sets were used. Above 40 EeV this corresponds to 
A^obs = 99 "AGASA+SUGAR" events, or to iVobs = 57 "AGASA only" events. Above 10 EeV comparison with an isotropic 
distribution of 1500 events was made, see text for more details. The last two likelihoods are for the auto-correlation Eq. I^J for 
9 < 20°, and the clustering within 2.5° up to multiplicity 10, respectively. The likelihoods are computed for n = 4 in Eq. @ 
which leads to reasonable discriminative power. 



# 


Us [Mpc 




structure 


observer 


-Bobs/G 




'-K10 


'-1 = 1 


^40 
'-t<10 


'-1=1 


^10 
'-K10 


rlO 
'-1=1 


'-B<20'> 


/-40 

'-n<10 
















AGASA only 




AGASA- 


hSUGAR 




AGASA only 


1 


2.4 X 10" 


-4 


yes 


1 


1.3 X 10"^ 


2.4 


0.070 


0.011 


0.37 


0.094 


0.12 


0.042 


0.57 


0.85 


2 


2.4 X 10" 


-4 


yes 


2 


8.2 X 10~^^ 


2.4 


0.43 


0.35 


0.52 


0.48 


0.16 


0.18 


0.52 


0.85 


3 


2.4 X 10" 


-4 


yes 


1 





2.6 


0.23 


0.15 


0.37 


0.39 


0.15 


0.15 


0.42 


0.73 


4 


2.4 X 10" 


-5 


yes 


1 





2.6 


0.25 


0.21 


0.33 


0.48 


0.11 


0.19 


0.30 


0.65 


5 


2.4 X 10" 


-5 


no 


1 





2.6 


0.36 


0.34 


0.45 


0.51 


0.13 


0.24 


0.65 


0.71 


6 


2.4 X 10" 


-5 


yes 


2 


8.2 X 10"^^ 


2.4 


0.49 


0.32 


0.79 


0.62 


0.17 


0.24 


0.56 


0.83 


7 


2.4 X 10" 


^4 


no 


1 





2.6 


0.35 


0.40 


0.42 


0.47 


0.12 


0.17 


0.53 


0.78 


8 


2.4 X 10" 


-6 


no 


1 





3.0 


0.36 


0.45 


0.19 


0.51 


0.10 


0.17 


0.24 


0.48 


9 


2.4 X 10" 


^4 


yes 


2 





2.6 


0.32 


0.31 


0.51 


0.49 


0.13 


0.20 


0.50 


0.75 


10 


2.4 X 10" 


-5 


yes 


2 





2.6 


0.24 


0.70 


0.32 


0.46 


0.10 


0.18 


0.46 


0.65 



of such peaks increases with the source density. For the 
same scenario 6, Fig.Elshows the distribution of UHECR 
deflection angles a with respect to the line of sight to the 
sources above various energy thresholds. This shows that 
deflection can be substantial even at the highest energies. 
Qualitatively this can be understood as follows: Neglect- 
ing energy loss processes, the r.m.s. deflection angle over 
a distance r in EGMF of r.m.s. strength B and coher- 
ence length Ic is 6{E,r) ~ (2r/c/9)i/2/rL ^j, where the 
Larmor radius of a particle of charge Ze and energy E is 
— E/{ZeB). In numbers this reads 




for r ^ Ic- Keeping in mind that sources are correlated 
with relatively strong fields and that according to Fig. [3 
~ 10% of the volume is filled with fields > 10 nG, we 
can see that deflections of order 20 degrees up to 10^" eV 
should therefore not be surprising. It is interesting to 
note in this context that, as can be seen from Figs. 12141 
the EGMF in our simulations tend to be more extended 
than the baryons and thus the distribution of sources if 
they follow the baryons. A significant part of the total 
deflection can therefore be contributed by fields that are 
not in the immediate environment of the sources. Fur- 
thermore, for the source density ~ 2.4 x 10~^Mpc~'^ 
of the typical scenario 6 of Tab. ^ the average distance 
to the closest source is ~ 15 Mpc, and thus close enough 
to make magnetized regions of several Mpc around the 



sources extend of order ten degrees on the sky, see also 
Fig. 01 upper panel. 

The analytical approximation for the time delay cor- 
responding to the deflection Eq. © is t{E, r) ~ 
re{E,r)^/4:, or 

These numbers are consistent with Fig. |S1 if one takes 
into account that the typical propagation distance drops 
from a few hundred Mpc at ^ 10^^ eV to ^ 30 Mpc above 
~ 10^° eV. Finally, Fig. [7| shows, for this same scenario, 
the predicted structure of arrival direction distributions. 

Sufficient data to definitely discriminate among the dif- 
ferent scenarios presented in Tab. HJ e.g. the presence 
of diffuse EGMF as well as UHECR source characteris- 
tics, will have to await the next-generation experimental 
facihties, such as the Pierre Auger ^ and EUSO |2^ 
projects. To demonstrate what can be achieved with 
these new experiments, we will work out the predictions 
for various observable statistical quantities for UHECRs 
both above 4 x 10^^ eV and lO^*^ eV for the Auger full sky 
observatory. For this we will assume an exposure about 
25 times larger than AGASA, which is reachable after 
about three years of running experiment. It corresponds 
to ~ 1500 events observed above 4 x 10^^ eV and ~ 30 
events above 10^" eV if a GZK cutoff is present, or ^ 200 
events above 10^° eV if no GZK cutoff is present. We will 
also sometimes consider ~ 5000 events above 4 x 10^® eV, 
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40. < E/EeV < 9000.: a= 63.B± 4S.a± 43.1 
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FIG. 6: The cumulative distribution of UHECR deflection 
angles with respect to the line of sight to the sources larger 
than a , averaged over 24 realizations for scenario 6 in Tab.Q] 
Shown are the distributions (middle, histogram) and 1-a vari- 
ations (upper and lower curves) above 4 x 10^^ eV (upper 
panel), above 10^" eV (middle panel), and above 2 x 10^° eV 
(lower panel). Also given on top of the figures average and 
variances of the distributions. 
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FIG. 7: Arrival direction distributions predicted within sce- 
nario 6 of Tab. 13 averaged over 34 realizations of 10 simu- 
lated trajectories above 10^® eV each. The filled red sphere 
represents the overall direction to the supergalactic centre. 
Upper panel: ~ 10^ events above 4 x 10^^ eV, convoluted 
with an angular resolution of 2.5°. Lower panel: ~ 10'' events 
above lO'^" eV, convoluted with an angular resolution of 1°. 



as is expected to be easily achieved by the EUSO exper- 
iment 23] . 



A. Signatures of structured versus homogeneous 
source distributions 

One expects that the nonuniformity in the source dis- 
tribution (structure) mostly influences the large scale 
multi-poles. Fig. |S1 shows that, for a given source den- 
sity, if magnetic deflection can be neglected, a homoge- 
neous distribution of sources predicts, not surprisingly, 
an angular power spectrum more isotropic than a struc- 
tured source distribution. Unfortunately, cosmic vari- 
ance is sufficiently large that one can definitely discrimi- 
nate the case of a structured distribution of sources only 
if the measured multi-poles happen to be several stan- 
dard deviations above the prediction of a homogeneous 
distribution. This can happen only when the source lu- 
minosities are almost uniform, which is rather unrealis- 
tic. Otherwise, the measured multi-poles are consistent 
with both structured and unstructured source distribu- 
tions, see Fig.|Sl In addition, within cosmic variance, the 
multi-poles depend insignificantly on the source density. 
In contrast, magnetic fields tend to significantly decrease 
their values, while at the same time reducing cosmic vari- 
ance, see Fig. O As a result, predictions of low-scale 
multi-poles above 4 x 10^^ eV for next generation exper- 
iments tend to deviate from isotropy more significantly 
than for unmagnetized sources. 

To demonstrate these tendencies more quantitatively, 
in Tab. UTI we provide the predictions for the sum of the 
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FIG. 8; Angular power spectrum C{1) as a function of 
multi-pole /, for the full sky Pierre Auger observatory as- 
suming A'obs = 5000 events observed above 4 x 10^^ eV. 
The upper panel assumes sources with identical properties, 
whereas the luminosity function Eq. Q was used in the lower 
panel. We show the realization averages (diamonds), statisti- 
cal (left) and total (right) error bars, respectively, predicted 
by the model. The red (thick) diamonds and red (thick, 
outer) error bars represent scenario 4 (unmagnetized struc- 
tured sources), whereas the black (thin) diamonds and black 
(thin, inner) error bars represent scenario 5 (unmagnetized 
homogeneously distributed sources). The source density is 
Us = 2.4 X 10~^ Mpc~^ in both cases. The straight line is 
the analytical prediction, C; ~ (47rA'obs)~^, for the average 
multi-poles for complete isotropy. 



first ten multi-poles for the scenarios from Tab.^ 



B. Source Density 

If magnetic field effects can be neglected, the small- 
scale auto-correlation function will depend strongly on 
the source density: Few sources imply strong auto- 
correlations with considerable cosmic variance as well 
as strong clustering, whereas many sources imply weak 



FIG. 9: Same as Fig. |H1 but comparing predictions of struc- 
tured sources with same densities and different magnetiza- 
tions. The red (thick) diamonds and red (thick, outer) error 
bars represent scenario 6 (magnetized sources), whereas the 
black (thin) diamonds and black (thin, inner) error bars rep- 
resent scenario 4 (unmagnetized sources). 



auto-correlation with comparatively small cosmic vari- 
ance. In the case of homogeneously distributed sources 
with identical properties this was indeed suggested as a 
measure of the source density ,19|]. Fig. ^| demonstrates 
that this is also true for structured sources. Furthermore, 
for a given source density the small scale auto-correlation 
is relatively insensitive to the structure of the source dis- 
tribution, as becomes evident in Fig. 1111 

However, if the sources are immersed in magnetic fields 
of order ^ 0.1/iG, as they can occur in galaxy clusters 
and filaments, the auto-correlation function becomes al- 
most insensitive to the source density and instead be- 
comes a probe of source magnetization, as discussed be- 
low in Sect. 3.D. This is because a structured EGMF 
of such strength diffuses cosmic rays up to 10^° eV over 
the whole region immersed in such fields, as can be seen 
from the rough estimate Eq. Therefore, the number 
of sources within such regions doesn't significantly affect 
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TABLE II: Predictions for the sum of the first ten multi-poles, 
^i-' above 4 X 10^^ eV for A^obs = 5000 observed events, 
with and without source luminosity variations according to 
Eq. lHJ. The two errors represent statistical and total (in- 
cluding cosmic variance) errors, as in the figures. The rows 
represent the scenarios of Tab.|I] 



2 6.87±0.39±5.9 2.21±0.16±0.58 

3 10.6±Q.52±17 4.48±Q.27±2.3 

4 10.8±Q.52±16 7.75±Q.43±8.5 

5 8.24±0.46±9.0 4.30±0.28±1.7 

6 4.11±0.25±2.7 2.27±0.17±0.69 

7 8.30±0.46±9.4 2.57±0.18±Q.55 

8 14.8±Q.75±16 14.9±0.73±16 

9 9.07±0.49±13 2.88±0.20±0.57 

10 11.4±0.58±12 8.38±0.48±5.1 
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FIG. 10: The auto-correlation function for A'^^bs ~ 200 events 
observed above 10^° eV for the full sky Pierre Auger ob- 
servatory for a bin size of A9 = 1°. Averages and er- 
rors are as in Fig. |H] Compared are unmagnetized struc- 
tured sources with different densities: The red (thick) dia- 
monds and red (thick, outer) error bars represent scenario 
4 {us = 2.4 X lO'^^Mpc"^), whereas the black (thin) dia- 
monds and black (thin, inner) error bars represent scenario 3 
{us = 2.4 X 10"" Mpc"^) in Tab. El 

the auto-correlations any more, which increases the un- 
certainty in the source density. 

C. Source Luminosity Function 

If the source luminosity function is parametrized by 
dus/dQ cx Q^" in an interval Qmin < Q < Qmax, then for 
a <^2 the most luminous sources will dominate, whereas 
for a ^ 2 the weakest sources will dominate. These two 
cases will therefore be approximated by scenarios with 
sources of identical properties and cosmic variance is ex- 
pected to be minimal. In contrast, a ~ 2 implies that all 
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FIG. 11: Same as Fig. 1101 but comparing unmagnetized 
sources with the same density Us = 2.4 x 10~^ Mpc~'^ with 
structured and unstructured distributions: The red (thick) di- 
amonds and red (thick, outer) error bars represent scenario 4 
(structured source distribution), whereas the black (thin) di- 
amonds and black (thin, inner) error bars represent scenario 
5 (homogeneous source distribution) in Tab.|I] 



luminosities contribute approximately equally to the ob- 
served flux. This is close to the case we chose in Eq. ^ 
which therefore tends to maximize cosmic variance. As 
can be seen from Fig. ^1 in this case the uncertainties 
of large scale multi-poles due to cosmic variance are in 
general larger than the ones due to the finite number of 
events observed, even for the relatively sparse data set 
currently available. For the small scale properties de- 
scribed by the auto-correlation function and cluster fre- 
quencies, for A^obs ^ 100, cosmic variance is in general 
smaller than or comparable to the statistical error, as 
can be seen from Figs. ^1 and ^2 and Figs. ^] and T\M 
below. 



D. Source Magnetization 

Scenarios 3 and 4 are somewhat disfavored by auto- 
correlation and clustering: Structured sources produce 
more clustering in the absence of magnetic fields, as can 
be seen, for example, in Fig. ^| While this effect is 
marginal in the current data set, it will develop into 
a strong discriminator for future experiments: Scenar- 
ios with no significant magnetic fields predict a stronger 
auto-correlation at small angles, independent of whether 
or not the sources are structured. This is because if 
sources are immersed in considerable magnetic fields, 
their images are smeared out, which also smears out the 
auto-correlation function over a few degrees. 

As shown in Fig. 1141 this effect is significant once 
^obs ^ 10'^ events are observed above 4 x 10^^ eV, as 
will be the case for the Pierre Auger project, and also if 
^obs ^ 100 events are observed above 10^° eV, as could 



luminosity function Eq. JTJ, constant luminosity, ^ 
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FIG. 12: The angular power spectrum C{1) as a function 
of multi-pole predicted for the AGASA+SUGAR exposure 
function (see text), for A'^nbs = 99 events observed above 40 
EeV, sampled from the simulated configurations of scenario 
2 in Tab. The diamonds indicate the realization averages, 
and the left and right error bars represent the statistical and 
total (including cosmic variance due to different realizations) 
error, respectively, see text for explanations. The histogram 
represents the AGASA+SUGAR data. The overall likelihood 
significance for n = A and / < 10 in Eq. ^ appears at the 
top of the figures. The upper panel takes into account cosmic 
variance due to variation of source properties assumed to be 
parameterized by Eq. The lower panel assumes that all 
sources have identical properties and thus cosmic variance is 
uniquely due to variation in source location. 



be achieved, for example, by the EUSO project. At these 
energies the strongest auto-correlation at small angles 
is predicted by scenario 4 in Tab. ^ where sources fol- 
low the large scale structure in the absence of magnetic 
fields. For a given source density, the auto-correlation 
predicted by homogeneously distributed sources in the 
absence of magnetic fields (scenario 5 in Tab. lies in 
between the one predicted by structured sources in the 
absence of magnetic fields (scenario 4 in Tab. P) and the 
one of structured sources with magnetic fields (scenario 



FIG. 13: The angular correlation function N{6) as a function 
of angular distance 6, predicted for a bin size of A6 = 1° 
for the AG AS A exposure function (see text), for A'obs = 57 
events observed above 40 EeV, sampled from the simulated 
configurations. Note that an isotropic distribution would cor- 
respond to N{6) = 1. Averages and errors are as in Fig. |H] 
The histogram represents the AGASA data. The two scenar- 
ios shown correspond to structured sources of equal density 
with (scenario 6, upper panel) and without (scenario 4, lower 
panel) magnetization. The overall likelihood significance for 
n = 4 and 6 < 10° in Eq. appears at the top of the figures. 
It is not significantly different for somewhat larger bin sizes 
Ae ~ 2°. 



6 in Tab.llll. 

As an example, the predictions for the auto-correlation 
function in the first degree bin are summarized in 
Tab. mil This also confirms the discussion in Sect. 3.B: 
For structured sources in the absence of magnetic fields, 
the degree-scale auto-correlation function decreases with 
increasing source density, whereas it is rather indepen- 
dent of the source density if the sources are immersed in 
magnetic fields. Finally we observe that the low-^ auto- 
correlation in case of magnetized sources hardly depends 
on their densities, as can be seen by comparing scenar- 
ios 2 and 6 in Tab. IIIII In contrast, for unmagnetized 
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FIG. 14: The auto-correlation function for A'obs ~ 1500 
events observed above 40 EeV for the full sky Pierre Auger 
observatory. Upper panel: Comparing magnetized sources 
with different densities. The red (thick) diamonds and red 
(thick, outer) error bars represent scenario 6 {us = 2.4 x 
10~^Mpc~'^), whereas the black (thin) diamonds and black 
(thin, inner) error bars represent scenario 2 {us — 2.4 x 
lO'^Mpc"^) in Tab. 13 Lower panel: Same as upper panel, 
but in the absence of magnetic fields, i.e. comparing scenario 
4 (red, thick set) with scenario 3 (black, thin set) in Tab.|I| 



sources the degree-scale auto-correlation is quite sensitive 
to the source density and always significantly higher than 
for magnetized discrete sources. This is clearly demon- 
strated by Fig. A consequence of this is that in the 
absence of magnetic fields the auto-correlation function 
and its fluctuations can indeed be used as a measure of 
the source density, as suggested in Ref. J^] : The strength 
of the auto-correlation at small scales as well as its cos- 
mic variance at all scales increase with decreasing source 
density. However, this effect can be almost completely 
erased by magnetic fields surrounding the sources. Such 
fields also tend to considerably reduce the effects of cos- 
mic variance, especially at low source densities, as Fig. ll4l 
demonstrates. This is also reflected by the fact that os- 
cillations predicted to occur in the actually measured 



TABLE III: Predictions for the auto-correlation function in 
the first degree bin, A'^(l°), for various threshold energies and 
number of observed events A^obs. The two errors represent 
statistical and total (including cosmic variance) errors, as in 
the figures. The rows represent the scenarios of Tab.Q] 





E > 10^° eV, 


S > 4 X 10^'^ eV, 


# 


Aobs = 200 


Aobs = 1500 


2 


45.4±12±17 


10.1±0.86±1.6 


3 


72.9±17±35 


16.3±1.6±4.3 


4 


98.0±21±53 


19.7±2.0±5.5 


5 


83.0±18±34 


15.9±1.5±3.4 


6 


44.8±12±19 


11.2±0.97±1.6 


7 


60.8±15±30 


13.1±1.3±2.2 


8 


191±33±88 


37.9±3.1±11.6 


9 


68.6±17±37 


14.4±1.4±3.6 


10 


100±21±40 


20.2±1.9±4.4 



auto-correlation function, i.e. predicted for a given re- 
alization of sources and magnetic fields, are considerably 
suppressed by magnetic fields. This is demonstrated by 
Fig.Uni 

Furthermore, Tab. Illll shows that the small-scale auto- 
correlation is rather independent of the observer position: 
The pairs of simulations 3, 9 and 4, 10 with unmagnetized 
sources differ only in the observer position, respectively, 
see Tab.m Their predictions for the average degree-scale 
auto-correlations indeed differ by amounts much smaller 
than even the statistical error. This suggests that the 
concrete realization of the local large scale structure is 
mostly irrelevant to tendencies presented in the current 
paper. 

As discussed in Sect. 3 A, magnetized sources also tend 
to suppress large scale multi-poles as well as their cosmic 
variance. However, since unmagnetized sources, in par- 
ticular in case of structured distributions, predict cosmic 
variances nearly as large as the averages of the multi- 
poles, small multi-poles cannot be used as a signature of 
strongly magnetized sources. 

E. Magnetic Fields surrounding the Observer 

If the observer is immersed in ~ O.l/xG fields, consid- 
erable \ow-l multi-poles are predicted. This possibility 
is already disfavored by current data as seen from 
Tab. n We demonstrate this by comparing the predic- 
tions of scenarios with a strongly or negligibly magne- 
tized observer for the large scale power spectrum above 
10^^ eV in Fig. 1161 This can qualitatively be under- 
stood as follows [23]: The presence of a region with a 
relatively strong magnetic field surrounding the observer 
preferentially shields off UHECRs from sources that are 
farther away from the observer, in particular those out- 
side of the magnetized region. In contrast, the flux from 
sources within the magnetized region is enhanced by the 
increased confinement time. Thus the observed flux is 
dominated by a few closer sources and appears more 



14 



N„ta( 100. < E/EeV < 1000.)= 200 




N„b,( 10. < E/EeV < 1000.)= 1500 



40 60 80 

theta (deg) 



0.0010 - 



0.0001 



-l — I — I — I — l- 



-l — I — I — I — l- 



J , , , L 



1 1* 
I 



I « 



J , , , L_ 



8 10 



100. < E/Eev < 1000.)= aoo 



100 




40 60 80 

theta (deg) 



FIG. 15: The auto-correlation function for the full sky Pierre 
Auger observatory, for A'obs ~ 200 events observed above 100 
EeV, predicted by one particular source realization for Us — 
2.4 X 10~^Mpc~^. Upper panel: scenario 6 in Tab. |l] i.e. 
with magnetic fields. Lower panel: scenario 4 in Tab. i.e. 
without magnetic fields. 



FIG. 16: The angular power spectrum obtained for the com- 
bined AGASA-f SUGAR exposure function, for A^'obs = 1500 
events observed above 10 EeV. The red (thick) diamonds and 
red (thick, outer) error bars represent scenario 6, whereas the 
black (thin) diamonds and black (thin, inner) error bars rep- 
resent scenario 1. The line key is as in Fig.|S] 
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anisotropic. 

Furthermore, scenarios in which the observer is im- 
mersed in an EGMF < 0.1^ G predict UHECR spectra 
with a pronounced GZK cutoff. In contrast, if the ob- 
server is strongly magnetized the cutoff is attenuated. 
One can understand this as a partial compensation be- 
tween energy attenuation lengths which decrease with 
increasing energy and magnetic diffusion lengths which 
increase with energy. As long as the source density is suf- 
ficiently high around the observer, the number of sources 
contributing to the flux decreases more slowly with in- 
creasing energy than in the absence of magnetic fields. 
For a strongly magnetized observer the observed spec- 
trum appears to be less modified with respect to the 
injection spectrum. This tendency is confirmed by our 
Monte Carlo simulations, see Fig. [T7I 

By comparing auto-correlations and clustering in sce- 
narios with and without EGMF we also confirmed that 
magnetic lensing j40| is insignificant for these highly 



FIG. 17: Predicted spectrum observable by AGASA for sce- 
nario 2 (weakly magnetized observer), as compared to the 
AGASA (dots) and HiRes-I 3 (stars) data. The histogram 
marks the average and the two lines above and below the 1- 
sigma fluctuations over the simulated realizations. The solid 
straight line marks the injection spectrum. For comparison 
shown (red, thin histogram) is the spectrum predicted by 
scenario 1 (strongly magnetized observer). Both cases cor- 
respond to the same source density Us — 2.4 x 10"* Mpc~"^. 

structured EGMF, even for cases where the observer is 
immersed in strong magnetic fields. 



IV. DISCUSSION AND CONCLUSIONS 

We studied the effects of non uniform UHECR source 
distributions and EGMF imprinted in the multi-pole mo- 
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ments, auto-correlation function, and cluster statistics of 
UHECRs with energies above 10^^ eV. We compared sce- 
narios where the sources are distributed homogeneously 
or according to the baryon density distribution obtained 
from a simulation of the large scale structure of the Uni- 
verse. We also compared the case in which the sources 
are all identical and in which their power and spectra 
of injected UHECRs are distributed according to func- 
tions that maximize cosmic variance. The changes in the 
results due to variation in the locations and properties 
of the sources were evaluated as cosmic variance. The 
source luminosity function was chosen close to the crit- 
ical case where all luminosities contribute comparably 
to the observed flux in a logarithmic distribution, which 
tends to maximize cosmic variance. In this case we found 
that cosmic variance of large scale multi-poles tends to 
be larger than the statistical fluctuations if the number 
of observed events is > 100. In all other cases, notably 
for the auto-correlations and cluster multiplicities, cos- 
mic variance is small compared to statistical fluctuations 
for iVobs ^100 events. 

The influence of EGMF is assessed by comparing the 
case of UHECR propagation when (a) EGMF are neg- 
ligible and (b) they are not and are modeled through 
a simulation of large scale structure formation and mag- 
netic field evolution. In the latter case the magnetic fields 
were normalized so as to reproduce published values of 
Faraday rotation measures for clusters of galaxies and fil- 
ametns ^^,"29^. We chose two cases for the location of the 
observer: in a relatively high field region with B ~ Cl/iG 
and in a negligible field region with B ~ 10^^^ G. 

In Ref. |3j it was found that in the case of struc- 
tured magnetic fields an observer immersed in ~ 0.1 fj,G 
fields was marginally favored whereas isotropy observed 
at 10^^ eV could not be explained by the local component 
alone. In Ref. 22] and in the present paper we took into 
account the contribution of sources beyond the large scale 
structure simulation box by making use of the periodic 
boundary conditions of the simulation box. As a result 
we found that the AGASA data is consistent with and 
even marginally favors scenarios with magnetic fields up 
to a few micro Gauss in galaxy clusters whose structure 
is obtained from the large scale structure simulation, and 
where the sources follow the baryon distribution. Con- 
trary to Ref. [3 this scenario requires that the observer 
is in a relatively low field region, B ^ 0.1 /iG. This also 
allows to explain isotropy observed at 10^^ eV reasonably 
well. That the fit quality there is significantly worse than 
at higher energies, as can be seen from Tab.^ may have 
to do with significant evolution of the characteristics of 
sources at large distances whose contribution increases 
with decreasing UHECR energy. The best fit scenario 
also predicts a pronounced GZK cutoff as well as consid- 
erable deflection of order 20° up to ~ 10^° eV. Thus, if 
this scenario is further confirmed by future experiments, 
charged particle astronomy may not be possible. 

Structured sources and magnetic fields up to a few mi- 
cro Gauss in galaxy clusters and the Earth immersed 



in relatively low EGMF B <^ 0.1 ^G in fact seem to 
be the most realistic scenario. Thus, already existing 
UHECR data allow probing the large scale structure dis- 
tribution of UHECR sources and magnetic field! We have 
demonstrated how future detectors such as the Pierre 
Auger and EUSO projects can further probe EGMF and 
UHECR source characteristics: Strongly magnetized ob- 
servers predict considerable large scale anisotropics be- 
tween 10^^ eV and 10^° eV which is already ruled out by 
current data on the percent level. Furthermore, strong 
fields surrounding the observer would predict a GZK 
cutoff that is less pronounced than for negligible fields. 
Low auto-correlations at degree scales imply magnetized 
sources quite independent of other source characteristics 
such as their density. The latter can only be estimated 
from the auto-correlations halfway reliably if magnetic 
fields have negligible impact on propagation. The multi- 
poles for I < 10 are relatively insensitive to the source 
density and on average tend to be smaller if sources are 
magnetized. At the same time, magnetized sources also 
tend to reduce cosmic variance of these multi-poles, and 
as a result their predictions above 4 x 10"'^^ eV for future 
experiments tend to deviate from isotropy more signif- 
icantly than for unmagnetized sources. Unfortunately, 
especially if source luminosities fluctuate considerably, it 
may be difficult to distinguish between structured and 
homogeneous source distributions even with next gener- 
ation experiments. 

For the required average local source number density 
and continuous power per source above 10^^ eV we find 
Us > 10^5 - 10-*Mpc"^ and Qs < lO-^^ergs^^ re- 
spectively, the latter within about one order of mag- 
nitude uncertainty to both sides. This corresponds to 
an average UHECR emissivity of (?uhecr — n-sQs ^ 
1.5 X lO'^^ergMpc"^ s""'^, with an uncertainty likely some- 
what smaller than for the above quantities, since it is 
fixed by the observed UHECR flux. Note that the uncer- 
tainty in Us is increased by the structured EGMF which 
tends to mix contributions from individual sources resid- 
ing in these fields. 

A simple study of deflection angles in the context of a 
constrained large scale structure simulation has recently 
been undertaken in Ref. [13]. They find maximal deflec- 
tion angles of a few degrees above 4 x 10^^ eV, consid- 
erably smaller than in our present study. This may be 
due to at least two reasons: First, these authors do not 
study structured sources which tend to be in regions of 
high density and magnetic fields which introduces a bias 
towards small deflection in the case of Ref. 17]. In fact, 
if we take the EGMF scenario corresponding to scenario 
2 in Tab. ^ but with a homogeneous source distribu- 
tion instead, we obtain average deflection angles of ~ 61° 
above 4 x 10^^ eV, ~ 33° above lO^" eV, and ~ 10° above 
2 X 10^° eV. This is smaller than the deflection angles 
obtained in scenarios where source positions and strong 
magnetic fields are correlated, see Fig. El but still consid- 
erably larger than values obtained in Ref. |0]. In fact, 
even if the magnetic field strength is reduced by a factor 
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10 in our simulations, the average deflection angle above 
4x 10^9 eV is still - 30°, only a factor ~ 2.2 smaller. This 
non-linear behavior of deflection with field normalization 
is mostly due to the strongly non-homogeneous character 
of the EGMF. Note that a field strength normalization 
reduction by a factor 10 corresponds to a left-shift of the 
x-axis m Fig. U] of two decades, such that the EGMF 
is virtually dynamically unimportant in all cells of the 
simulation. 

The remaining discrepancy in typical UHECR deflec- 
tion between our present work and Ref. for cases 
meaningful to compare, is most likely due to the dif- 
ferent EGMF models used by the different authors. A 
detailed consideration of these differences is well beyond 
the scope of this paper, although we plan to investigate 
this issue further in a forthcoming paper. So little is 
known about large scale magnetic fields and their evolu- 



tion. This makes signatures for magnetized sources, as 
discussed in the present paper, even more important. 
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